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I .  INTRODUCTION 

Probahly  the  first  contribution  of  missile  stability  analysis  was 
the  placement  of  feathers  on  the  arrow.  This  created  a  restoring  moment 
and  made  the  arrow  a  statically  stable  missile  whose  angular  motion 
is  approximated  by  a  sine  wave.  Early  bullet,  designers  soon  realized 
that  feathers  or  fins  are  most  inconvenient  for  gun  launch  and 
found  it  necessary  to  impart  a  high  spin  rate  to  the  bullet  or  shell 
by  means  of  rifled  tubes.  If  sufficient  spin  is  given  the  projectile, 
it  is  gyroscopically  stable  and  its  angular  motion  is  approximated  by 
the  sum  of  two  sine  waves.  From  a  dynamics  point  of  view,  the  motion 
of  a  symmetric  statically  stable  missile  is  a  special  case  of  the 
motion  of  a  gyroscopically  stable  missile. 

Until  World  War  I,  the  primary  concern  of  the  designer  was  static 
or  gyroscopic  stability  and  the  most  important  moments  were  the  static 
moments,  which  were  assumed  to  be  linear  in  the  angle  of  attack  or 
sideslip: 
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If  is  not  zero,  the  aerodynamic  moment  of  Eq.  (3)  is  essentially 
a 

asymmetric  and  the  missile's  angular  motion  lias  the  complexity  of 

aircraft,  motion.  If  both  Cw  and  £  are  zero,  the  aerodynamic  moment 

M  <L 

a 

is  symmetric  and  the  resulting  angular  motion  is  that  of  a  body  of 

revolution.  For  the  intermediate  case  of  zero  C„  and  nonzero  e  , 

M  a 


; 


i 


7 


the  missile  has  a  slight  aerodynamic  asymmetry  and  its  moment  is  that 
of  a  body  of  revolution  with  respect  to  the  complex  aerodynamic  trim 
angle, 

Since  World  War  II,  missile  designers  have  encountered  a  number 
of  surprises  with  regard  to  the  dynamic  stability  of  their  designs. 

The  dynamic  stability  is  influenced  by  additional  moment  terms  and 
determines  the  growth  or  decay  of  the  oscillatory  angular  motion. 

In  this  paper,  we  will  give  a  survey  of  dynamic  instabilities  that 
have  been  observed  as  well  as  some  that  are  possible  but  not  yet 
observed.  For  most  of  this  paper,  only  symmetric  missiles  or  missiles 
with  slight  asymmetries  will  be  considered.  A  brief  discussion  of  an 
almost  symmetric  missile  ( |C^  |<<|C^  |)  as  well  as  the  effect  of  moving 

payloads  will,  however,  also  §e  given.  Equal  transverse  moments  of 
inertia  will  be  assumed  throughout  the  report. 


II.  NONSPINNING  SYMMETRIC  MISSILES 


In  addition  to  the  static  moment  coefi'icients,  the  linear  motion 
of  a  nonspinning  missile  is  affected  by  four  damping  moment  coefficients 


For  the  oscillatory  motion,  the  angular  derivatives  are  related. 

q  =  a  ,  r  =  -3  (41 

The  symmetry  assumption,  then,  allows  the  following  simple  expression 
for  the  complete  linear  transverse  moment: 
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The  angular  motion  for  a  statically  stable  missile  has  the  very  simple 
form 
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The  motion  is  a  damped 
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ellipse  with  semi -major 
It  is  important  to  note 


axis  K  +  K  and  semi- 
1  2 

that  A  =  A2  implies 


that  the  eccentricity  of  the  elliptical  motion  is  maintained.  Thus, 
initially  planar  motion  remains  planar  and  initially  circular  motion 
remains  circular. 


According  to  the  definition  of  A^ ,  the  motion  will  be  unstable  if 
H  is  negative.  The  primary  component  of  H  is  +  CM  and,  thus,  this 


M. 


instability  occurs  when  C 


M 


+  C„  is  positive.  In  other  words,  the 
a 


aerodynamic  moments  act  to  increase  the  angular  rates.  This  unexpected 
behavior  has  been  observed  at  hypersonic  speeds  for  slowly  spinning  re¬ 
entry  shapes.  This  dynamic  instability  could  be  caused  by  the  entropy 
gradient  induced  by  the  bow  shock  and  reinforced  by  ablation1.  At 
transonic  speeds,  unstable  damping  has  been  observed  by  MacAl lister2 
and  this  has  been  explained  by  nose-induced  separation-3. 


An  important  feature  of  MacAl lister's  ballistic  range  measurements 
was  that  the  initial  almost-planar  motion  quickly  became  an  oval 
almost-circular  limit  motion.  The  theoretical  explanation  of  this 
limit  motion  requires  the  introduction  of  a  rather  strange  damping 
moment.  If  we  rewrite  the  linear  damping  mompnt  terms  of  Eq .  (5)  using 
the  polar  form  of  the  angle  of  attack,  £  =  6e10  ,  we  have 


1.  L.E.  Ericsson ,  "Unsteady  Aerodynamics  of  an  Ablating  Flared  Body 
of  Revolution  Including  Effect  of  Entropy  Gradient ,  "  AIAA  Journal 
6j  December  1968 ,  pp.  2395-2401 . 

2.  L.C.  MacAllister,  "Some  Instability  Problems  with  Re-entry  Shapes ,  " 
Ballistic  Research  Laboratories  Memorandum  Report  No.  1224, 

August  1959,  AD  377344. 

3.  L.E.  Ericsson  and  J.P.  Reding ,  "Dynamic  Stability  Problems  Associ¬ 
ated  with  Flare  Stabilizers  and  Flap  Controls ,  "  Journal  of  Space¬ 
craft  and  Rockets  7_,  February  1970,  pp.  132-137. 
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CM  ’  -  4  tCM.  *  CM  5  tS'*  le's><>19  (») 

d  a  q 

Thus  for  linear  theory,  the  damping  moment  in  the  plane  of  the  total 
angle  of  attack  is  proportional  to  6 '  ,  the  radial  rate  of  change  of 
this  angle,  the  damping  moment  perpendicular  to  this  plane  is  pro¬ 
portional  to  6,  the  circumferential  rate  of  change  of  this  angle, 
and  the  proportionality  factors  are  equal. 

The  simplest  nonlinear  extension  of  this  damping  moment  expression 
is 


CM  =  -  i  (dQ  +  d262)  («'  +  i0/6)ei0  (9) 

d 

A  more  general  expression  that  retains  the  in-plane*  and  out-of-plane 
damping  equality  is  obtained  by  making  a  function  of  62  .  Even  this 
assumption  is  not  sufficient  to  generate  the  circular  limit  motion 
observed  by  MacAJ lister. 

A  successful  approach4  is  to  drop  the  equality  of  in-plane  and 
out-of-plane  damping: 


CM  -  -  i|d0  (&'  +  ie'6)  ♦  d262  [Cl  +  a) o'  +  ie'61  jei6 

=  -i  [ CdQ  +  d  62)  +  d2  a6$'£]  (10) 

For  constant  d2  and  a,  Eq.  (10)  introduces  two  cubic  damping  terms5. 
This  nonlinear  moment  expression  can  be  used  in  the  usual  quasilinear 


*An  "in-plane  moment "  means  a  moment  producing  a  rotation  in  the  plane 
of  the  total  angle  of  attack;  the  in-plane  moment  vector  is  thus 
normal  to  this  plane .  Similar  remarks  apply  to  the  out-of-plane 
moment. 

4.  C.H.  Murphy,  "Circular  Pitching  and  Yawing  Motion  of  Nose  Cone 

Configurations,  "  Ballistic  Missiles  and  Space  Technology  II, 
Pergamon  Press,  New  York,  1961,  pp.  328-336.  ( See  also  Ballistic 

Research  Laboratories  Report  No.  1071,  March  1959,  AD  216341.) 

5.  C.H.  Murphy,  "Slender  Body  Estimates  for  Too  Cubic  Aerodynamic 
Damping  Moments,  "  AIM  Journal  £,  March  1966,  pp.  536-537 . 
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analysis6-9.  According  to  this  analysis,  the  nonlinear  solution  can  be 
approximated  by  a  solution  of  the  form  of  Eqs.  (6-7)  in  which  the  X!s 

become  functions  of  the  K'.s,  9 

3 

X  =  -  [H  +  11  (K2  +  a  K2)  j  /  2  (11) 

10  2  12 


X  =  -  (H  +  H  (K2  ♦  a  K2)  ]/ 2  (12) 

2  0  2  2  1 

The  behavior  of  a  nonlinear  solution  can  be  described  by  trajectories 
in  a  K2  vs  K2  amplitude  plane.  Since  Eq .  (6)  for  X_.  =  0  generates 

ellipses,  each  point  in  the  amplitude  plane  identifies  an  elliptical 
motion  and  the  trajectory  through  that  point  describes  how  this 
elliptical  motion  changes  under  the  influence  of  nonlinear  damping. 
Points  on  the  amplitude  plane  axes  represent  circular  motions  and  the 
line  K2  =  K2  is  the  locus  of  planar  motions. 

If  11^  and  H  a^e  opposite  in  sign,  a  circular  singularity  exists 
with  amplitude  6 ^  -  [  -H  /H  ]*5  .  The  amplitude  plane  for  this  case 

and  for  equal  in-plane  and  out-of-planc  damping  (a  =  0)  is  given  in 
fig.  1.  The  circular  limit  motions  are  unstable  hut  there  is  a  stable 
planar  limit  motion  with  amplitude  2<$c .  It  can  be  shown  that  for 

a  <  1  the  circular  motions  are  unstable  but  for  a  >  1  they  are  stable. 


6.  C.H.  Murphy,  "The  Prediction  of  Nonlinear  Pitching  and  Yawing 
Motion  of  Symmetric  Missiles,  "  Journal  of  the  Aeronautical 
Sciences  24,  July  1957,  pp.  473-479.  (See  also  Ballistic  research 
Laboratories  Report  No.  99 5,  October  1956,  AD  122221.) 

7.  C.H.  Murphy ,  "Quasi-linear  Analysis  of  the.  Nonlinear  Motion  of  a 
Nonspinning  Symmetric  Missile ,  "  Journal  of  Applied  Mathematics  and 
Physics  ( 7AMP)  14 ,  No.  5,  September  25 ,  1963,  pp.  630-643. 

8.  W.H.  Haseltine ,  "Existence  Theorems  for  Nonlinear  Ballistics ,  " 

J,  Soa .  Indust,  Appl .  Math .  II,  September  1963 ,  pp.  553-563. 

9.  C.H.  Murphy ,  "Angular  Motion  of  a  Re-Entering  Symmetric  Missile , " 
AIAA  Journal  3,  July  1965 ,  pp.  1275-1282.  (See  also  Ballistic 
Research  Laboratories  Report  No.  1114,  August  1960 ,  AD  247271 ,  and 
Ballistic  Research  Laboratories  Memorandum  Report  No.  1358 , 

June  1961 ,  AD  266513.) 
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2  gives  the  amplitude  plane  for  a  =  2.  Numerical  integrations 
of  the  complete  equations  of  motion  verify  that  stable  circular  limit 
motions  do  exist  for  a  >1.  Recently,  several  authors10-13  have  made 
a  variety  of  wind  tunnel  measurements  of  out -of -plane  damping  and  have 
shown  that  it  can  be  quite  different  from  in-plane  damping  for  cones 
at  supersonic  and  hypersonic  speeds.  Since  1  +  a  is  the  ratio  of  the 
planar  damping  to  the  circular  damping,  we  see  that  this  ratio  must 
exceed  two  before  stable  circular  motion  can  exist. 

One  common  feature  of  Figs.  1-2  is  that  large-amplitude  planar 
motions  decay  and  small -amplitude  motions  grow.  Although  these  motions 
go  to  different  limit  motions,  the  final  motions  are  bounded.  Thus,  we 
could  expect  that  for  all  cases  when  planar  damping-in-pitch  wind 
tunnel  measurements  show  similar  behavior,  bounded  flight  motions  would 

occur.  Fig.  S,  for  H  and  H  both  negative  and  a  <  -1,  shows  that  all 

0  2 

motions  go  to  large  spiral  motions  although  planar  motions  tend  to  the 
planar  singular  motion  with  amplitude  2[(1  +  a)H2/HQ]-^  This  shows 

that  intuitive  arguments  should  be  applied  with  care  to  nonlinear 
systems14 . 


10.  M.  Tobak ,  L.B.  Schiff  and  V.L.  Peterson,  "Aerodynamics  of  Bodies 
of  Revolution  in  Coning  Motion ,  "  AIAA  Journal  7_,  January  1969 , 
pp.  95-99 . 

11.  L.B .  Sahiff  and  M.  Tobak,  "Results  from  a  New  Wind-Tunnel 
Apparatus  for  Studying  Coning  and  Spinning  Motions  of  Bodies  of 
Revolution,  "  AIAA  Journal  8,  November  1970,  pp.  1953-1958. 

12.  G.W.  Stone,  E.L.  Clark,  Jr.,  and  G.E.  Burt,  "An  Investigation 
of  Nonsyrmetric  Aerodynamic  Vamping  Moments,  "  AIAA  Paper  72-29, 
San  Diego,  Califoimia,  1972. 

13.  0.  Walchner  and  F.M.  Sawyer,  "'In-Plane'  and  'Out-of -Plane' 
Stability  Derivatives  of  Slender  Cones  at  Mach  14,"  Aerospace 
Research  Laboratories,  Wright-Patterson  Air  Force  Base ,  Ohio, 
Report  ARL  73-0090,  July  1973. 

14.  C.H,  Murphy,  "An  Erroneous  Concept  Concerning  Nonlinear  Aero¬ 
dynamic  Damping,  "  AIAA  Journal  1^  June  1963,  pp.  1418-1419. 
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III.  LINEAR  MAGNUS  MOMENT 
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For  a  statically  unstable  missile  such  as  a  shell  or  bullet,  high 
rates  of  spin  are  required  for  stability  and  a  Magnus  side  moment  must 
be  added  to  the  total  aerodynamic  moment.  Since  statically  stable 
missiles  are  usually  spun  to  reduce  the  effect  of  manufacturing  asym¬ 
metries,  this  Magnus  moment  should  be  considered  for  finned  missiles  as 
well  as  bodies  of  revolution. 

Before  doing  so,  wc  must  make  a  decision  on  the  appropriate  co¬ 
ordinate  system.  In  order  to  avoid  the  algebraic  complexities  of  a 
spinning  missile-fixed  XYZ  axis  system,  we  will  use  nonspinning  aero- 
ballistic  coordinates  XYZ.  The  X-axis  pitches  and  yaws  with  the 
missile.  The  Y-axis  is  selected  to  lie  in  the  horizontal  plane  ini¬ 
tially  but  it  is  very  important  to  note  that  it  does  not  remain  there. 
If  this  axis  were  required  to  remain  in  the  horizontal  plane,  the  spin 
of  the  system  would  be  nonzero  and  these  axes  would  be  called  fixed- 
plane  axes.  Fixed-plane  axes  are  useful  when  horizontal  or  vertical 
forces  such  as  gravity  are  present15.  Indeed,  it  has  been  shown  that 
an  ascending  or  descending  missile  with  a  constant  horizontal  control 
moment  could  have  an  instability  due  to  the  terms  involving  fixed-plane 
coordinates  spin16-17. 

In  aeroballistic  coordinates,  the  complex  angle  of  attack  is 
written  as  \  and  the  transverse  aerodynamic  moment  coefficients  become 


(13) 


15.  C.H.  Murphy,  "Gravity-Induced  Angular  Motion  of  a  Spinning 
Missile,  "  Journal  of  Spacecraft  and  Kocke ts  8j  August  1971, 
pp.  824-828.  (See  also  Ballistic  Research  Laboratories  Report 
No.  1546,  July  1971,  AD  730641.) 


16.  K 


17. 


si.H.  Lloyd  and  D.P.  Brown ,  "Instability  of  Spinning  Projectiles 
During  Terminal  Guidance,  "  Jouma l  of  Guidance  and  Control  2, 
January -February  1979,  pp.  65-70. 

C.H.  Murphy ,  "Instability  of  Controlled  Projectiles  in  Ascending 
or  Descending  Flight, "  AIAA  Paper  79-1669,  August  1979.  (See 
also  USA  ARFADCQM  Ballistic  Research  Laboratory  Memorandum  Report 
No.  02915,  April  1979,  AD  A072808.) 


The  first  term  on  the  right  side  of  Eq.  (13)  represents  the  Magnus 
moment,  which  is  proportional  to  the  spin  and  the  angle  of  attack.  A 
positive  Magnus  moment  coefficient  represents  a  moment  that  acts  to 
rotate  the  nose  of  the  missile  around  the  velocity  vector  in  the  direc 
tion  of  spin. 


The  presence  of  spin  and  Magnus  moment  does  not  change  the  form  of 
the  epicyclic  solution  of  Eq.  (6)  but  does  give  more  complicated  damping 
rates  and  frequencies. 


(14) 


[P  ±  /P2  -  4M] / l 


(15) 


Dynamic  stability  requires  positive  A j s  and  can  be  stated  simply  in 


terms  of  two  stability  factors, 


s  and  s,,  where 
g  d’ 


s  =  P2/4M 
g 


(16) 


sd  =  2T/H 


(17) 


The  gyroscopic  stability  factor,  s  ,  is  essentially  the  ratio  of 

g 

squared  gyroscopic  spin  to  static  moment  coefficient.  Periodic  motion 
occurs  when  this  stability  factor  is  greater  than  unity.  The  dynamic 


stability  factor,  s  , . ,  is  essentially  the  ratio  of  the  Magnus  moment 
coefficient  to  the  sum  of  the  damping  moment  coefficients.  For  dynamic 
stability18 


4M 

P2 


<  C2  -  Sd)5d 


(18) 


18.  C.H.  Murphy,  " Criteria  for  the  Generalized  Dynamic  Stability  of  a 
Rolling  Symmetric  Missile ,  "  Journal  of  the  Aeronautical  Sciences 
24j  October  1957 ,  pp.  773-774. 
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Fig.  4  summarizes  all  the  implications  of  this  relation.  Note 
that  if  s,  is  outside  the  interval  (0,  2),  a  statically  unstable 
missile  cannot  be  dynamically  stabilized  by  spin  and  a  statically 
stable  missile  can  be  made  dynamically  unstable  by  sufficiently  high 
rates  of  spin.  Very  simple  cone  cylinders  and  finned  cone  cylinders 
have  been  shown  to  have  s^  values  outside  this  region19.  The  linear 

Magnus  moment  is  largest  for  long  projectiles  with  boattails  at  tran¬ 
sonic  speeds.  Indeed,  one  155mm  developmental  shell  when  fired  at 
Mach  numbers  between  0.92  and  0.96  experienced  a  number  of  3  km  shorts 
due  to  Magnus  instability. 


IV.  NONLINEAR  MAGNUS  MOMENT 

In  1951,  the  Navy  was  faced  with  a  very  strange  problem20-21 
during  the  development  of  the  12.75"  antisubmarine  ship-launched  spin¬ 
ning  finned  rocket,  the  Weapon  A.  When  fired  to  the  port  side  of  a 
high-speed  destroyer,  it  performed  well.  When  fired  to  the  starboard 
side,  however,  its  angular  motion  grew  to  a  very  large  amplitude  coning 
motion  and  its  performance  was  completely  unsatisfactory. 

This  dependence  of  missile  stability  on  launch  conditions  is  a 
characteristic  of  nonlinear  differential  equations  and  the  cause,  in 
the  case  of  Weapon  A,  was  a  strongly  nonlinear  Magnus  moment.  The  be¬ 
havior  can  be  easily  predicted  by  the  quasil inear  theory  for  a  simple 
cubic  Magnus  moment,  i.e.  a  quadradic  Magnus  moment  coefficient: 


(19) 


For  this  moment,  the  quasil inear  exponential  damping  functions  which 
determine  trajectories  in  the  amplitude  plane  become 


19.  C.H.  Murphy ,  "Effect  of  Roll  on  Dynamic  Instability  of  Symmetric 
Missiles, "  Journal  of  the  Aeronautical  Sciences  21,  September 
1954,  pp.  643-644. 

20.  I.E.  Highberg ,  " Suggested  Mechanism  for  the  Instability  of  Weapon 
A,"  NOTS  TN-5036-94,  19  July  1951. 

21.  W.R.  Haseltine,  "Instability  of  Weapon  A  in  Cross  Winds,"  NAVORD 
Report  2057,  September  1953,  AD  0  2349  2. 
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(20) 


X1  = 


10 


+  f  b  c. 


CK2  +  2  Kp 


A  =  A 


20 


V  b  c. 


(Kj  +  2 


(21) 


A  possible  amplitude  plane  for  these  Ajs  is  given  in  Fig.  5.  According 

to  Eq.  (15),  the  frequencies  for  a  finned  projectile  (M  <  0)  are  oppo- 
site  in  sign  and  thus  its  complex  angular  motion  is  described  by  the 
sum  of  oppositely  rotating  two-dimensional  vectors.  Gravity  tip-off 
plus  the  crosswind  produced  by  launch  to  the  port  of  ships  moving  at 
thirty  knots  produces  10°  amplitude  clockwise  angular  motion.  Launch 
to  starboard  produces  10°  amplitude  counterclockwise  motion.  If  8^  =  5° 

the  motion  associated  with  port  launch  lies  to  the  left  of  the  dashed 
curve  (the  separatrix)  and  will  damp  to  a  small -amplitude  coning  motion 
Starboard  launch  is  to  the  right  of  the  separatrix  and  large -amplitude 
motion  is  successfully  predicted  by  the  theory. 


V.  AERODYNAMIC  TRIM  (E  t  0) 

a 


The  aerodynamic  moment  of  a  symmetric  missile  will  not  be  zero  at 

zero  angle  of  attack  if  its  body  or  fins  are  slightly  deformed  or  its 

center  of  mass  is  not  on  its  axis  of  symmetry.  This  trim  moment 

causes  the  £  in  Eq.  (3).  In  aeroballistic  coordinates,  the  resulting 
3- 

trim  angle  rotates  with  the  missile  and  has  an  amplitude  that  depends 
on  the  spin  rate22. 


+ 


+  E  k 
a  3 


i<j> 

e 


(22) 


k 


3 


-  <f/  +  i  A  ) 
yl  1  / 


(23) 


For  most  spins,  A  can  be  neglected.  Near 
important  and  determines  the  maximum  value 


resonance, 

of  |k  )  . 

3 


=  ‘I’j,  it  is 
A  plot  of  jk 3  | 


22.  J.D.  Niaolaidea ,  "On  the  Free  Flight  Motion  of  Missiles  Having 
Slight  Configurational  Asymmetries t  "  Ballistic  Research  Labora¬ 
tories  Report  No.  858j  June  1953t  AD  26405.  (See  also  IAS 
Preprint  395}  January  1953.) 
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versus  4/  is  given  in  Fig.  6  and  shows  the  need  for  missile  designers 
to  avoid  resonance. 


VI.  INDUCED  ROLL  MOMENT  (£  t  0) 

cl 

In  general,  roll  motions  are  such  that  a  missile's  roll  rate  will 
vary  through  resonant  spin  rate  and  the  missile  will  attain  only  a 
fraction  of  its  maximum  trim  value23'24.  If  the  roll  moment  at  angle 
of  attack  is  a  function  of  0,  the  roll  orientation  of  the  plane  of  the 
angle  of  attack25,  the  roll  moment  coefficient  can  be  written  in  the 
form 


C;  =  c£  +  *'  C  +6  C  (0,  6) 

0  'p  6 


(24) 


0  and  <5  can  be  computed  from  Eq .  (22). 


iO  i  (<£  1  -$)  i 

tie1  =  K  e  +  K2  e 


^ak3 


(25) 


For  most  motions,  0  varies  rapidly  and  the  rolling  motion  is  unaffected 
by  C  .  For  pure  trim  motion  (K,  =  K2)  or  two-mode  motion  near 
6 

resonance  (K  =  0,  <j>y _  =  4/ )  ,  0  is  constant  and  C  can  have  an  impor- 

di  1  )L 

O 

tant  effect.  Near  resonance,  6  can  be  quite  large  and  the  resulting 
induced  roll  moment  (i.e.,  the  third  term  in  Eq .  (24))  can  force  the 
roll  rate  to  stay  near  resonance.  This  phenomenon  of  "roll  lock-in" 
has  been  observed  in  flight  as  well  as  in  a  number  of  computer  simula¬ 
tions  . 


23.  R.J.  Tolosko,  "Amplification  Due  to  Body  Trim  Plane  Rotation , " 
AIAA  Paper  71-48j  January  1971. 

24.  C.H.  Murphy ,  "Response  of  an  Asymmetric  Missile  to  Spin  Varying 
through  Resonance,"  AIAA  Journal  9 ,  November  1971 ,  pp.  2197-2201. 
(See  also  Ballistic  Research  Laboratories  Report  No.  1242,  July 
1971 ,  AD  729772.) 

25.  J.D.  Ni'iolaides ,  "Two  Non-linear  Problems  in  the  Flight  Dynami-cs 
of  Modern  Ballistic  Missiles IAS  Report  59-17,  January  1959. 
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The  induced  roll  moment  can  be  caused  by  aerodynamic  asymmetry 
present  in  a  "symmetric"  four-finned  missile  or  by  mass  asymmetries  in 
an  aerodynamical ly  symmetric  missile.  For  example,  the  lift  force  on  a 
symmetric  missile  at  angle  of  attack  acts  on  the  center  of  mass  of  the 
missile.  If  the  center  of  mass  is  not  on  the  missile's  axis  of  sym¬ 
metry,  it  will  produce  a  roll  moment  of  the  form  of  Eq.  (24).  It 
should  be  emphasized  that  the  occurrence  of  "roll  lock-in"  depends  on 
the  details  of  the  pitching  motion  and  its  coupling  to  the  rolling 
motion  through  Eq.  (24)  and  can  only  be  determined  by  numerical  inte- 
grations26-  . 


VII.  INDUCED  SIDE  FOMENT  (£  t  0) 

Although  flight  failures  have  been  explained  by  the  occurrence  of 
resonance  through  "roll  lock-in,"  in  some  cases  angles  of  attack  have 
been  observed  much  larger  than  those  predicted  by  Eq.  (23) .  In  1959, 
Nicolaides25  developed  his  "catastrophic  yaw"  theory  by  the  introduction 
of  induced  side  moments.  The  existence  and  the  effect  of  these  moments 
have  been  discussed  hy  other  authors29-31.  Nicolaides'  induced  side 
moment  term  can  be  included  in  the  aerodynamic  moment  expression  of 
Eq.  (13)  by  adding 


26 .  L.  Glover,  "Effects  on  Roll  Rate  of  Maos  and  Aerodynamic  Asym¬ 
metries  for  Ballistic;  Re-entry  Bodies ,  "  Journal  of  Spacecraft 
and  Rocket 3  2,  Maroh-April  1965,  pp,  220-225. 

27.  D.A.  Price,  Jr.,  "Sources,  Mechanisms ,  and  Control  of  Roll 
Resonance  Phenomena  for  Sounding  Rockets,  "  Journal  of  Spacecraft 
and  Rockets  4,  November  1967,  pp.  1516-1525. 

28.  D.A.  Price ,  Jr.,  and  L.E.  Ericsson,  "A  New  Treatment  of  Roll- 
Pitch  Coupling  for  Ballistic  Re-Entry  Vehicles ,"  AIAA  Journal  £, 
September  1970,  pp.  1608-1615 . 

29.  F.J.  Regan ,  V.L.  Shermerhom  and  M.E.  Fatusi,  "Re1 '-l -Induced  Force 
and  Moments  Measurements  of  the  M823  Research  Store,  "  NOLTR 
68-195,  November  1968. 

30.  T.A.  Clare,  "Resonance  Ins'tability  for  Finned  Configurations 
Having  Nonlinear  Aerodynamic  Properties ,  "  Journal  of  Spacecraft 
and  Rockets  8,  March  1971 ,  pp.  278-283. 

31.  T.R,  Pepitone  and  I.D.  Jacobson ,  "Resonant  Behavior  of  a  Sym¬ 
metric  Missile  Having  Roll  Orientation-Dependent  Aerodynamics ,  " 
Jour’nal  of  Guidance  and  Control  tj  September-October  1978 , 

pp.  335-339, 
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[CSM  (8,  6)]  i 

a 


to  the  right  side  of  Eq.  (13). 


For  flight  conditions  for  which  0  is  constant,  the  primary  effect 
of  this  term  is  to  change  X  to 


H  <f/  -  P  T  + 


*1  *  - 


2  <t>i  -  P 


b  C,. 


(26) 


The  presence  of  the  induced  side  moment  coefficient  in  Eq.  (26)  intro¬ 
duces  the  possibility  of  a  very  small  X  ,  which  can  cause  a  very  large 

resonance  value  of  in  Eq.  (23).  An  even  worse  possibility  is  a 

large  positive  value  of  X^,  which  would  cause  an  exponential  growth  of 

k3  .  This  possibility  is  the  "catastrophic  yaw"  of  Nicolaides  and  may 

be  the  cause  of  some  spectacular  flight  failures. 


VIII.  NONLINEAR  AERODYNAMIC  MOMENT  (?  t  0) 

3. 

Tne  combination  of  nonlinear  aerodynamic  moments  with  a  trim 
moment  can  give  rise  to  a  rich  variety  of  limit  motions.  Many  of  these 
motions  have  been  produced  by  computer  simulations  but  as  yet  none  have 
been  observed  in  flight.  In  this  section,  we  will  briefly  consider  the 
effect  of  (a)  a  cubic  static  moment,  and  (b)  two  cubic  damping  moments. 

First,  the  static  moment  coefficient  is  assumed  to  have  the  form: 

Cm  +  1  Cn  =  "  1  £>0  +  C2  ^  C  -  C0  *a]  C27) 


For  pure  trim  motion,  this  assumption  replaces  the  linear  Eq .  (23)  by 
a  cubic  equation  for  k3  .  Near  resonance,  three  values  of  k3  can  be 

computed  but  only  two  correspond  to  stable  motion.  Much  more  interest¬ 
ing  limit  motions  have  been  found  which  are  generalized  subharmonic 
motions.  For  these  motions,  certain  constant  values  of  K  ,  K? ,  kg  have 

been  predicted  and  produced  by  computer  integrations  for  spins  far  from 


19 


resonance  as  well  as  quite  near  to  resonance32  33. 

A  second  set  of  limit  motions  can  be  constructed  for  the  damping 
moment  expansion  of  Eq.  (10): 


C 

m 


+  i  C 

n 


i  C 


M  «  -  V 
a 


+  C. 


(28) 


A  number  of  one-,  two-,  and  three-mode  limit,  motions  have  been  pre¬ 
dicted  by  the  quasil inear  theory  and  computed  by  numerical  integration 
of  the  equations  of  motion34. 


IX.  MOVING  INTERNAL  PARTS 

In  1955,  an  eight-inch  shell  showed  a  strange  spin  decay  coupled 
with  a  significant  range  loss35  (Fig.  7) .  The  range  loss  was  due  to  a 
growth  in  the  high-frequency  component  of  the  pitching  motion  (A1  >  0) . 

Recent  ballistic  range  tests  of  a  developmental  20mm  projectile  with 
the  M505  fuse  have  shown  a  growth  of  the  high-frequency  mode  and  un¬ 
explained  spin  decay36.  In  both  cases,  the  projectiles  carried  a 
component  that  could  move  a  small  amount  during  flight. 


32.  C.H.  Murphy,  "Generalized.  Subharmcnic  Response  of  a  Missile  with 
Slight  Configurational  Asytrrne tries ,  "  AIAA  Paper  72-972,  September 
1972.  (See  also  Ballistic  Research  Laboratories  Report  No.  1591, 
June  1972,  AD  749787.) 

33.  C.H.  Murphy,  B.A.  Nodes  and  J.W.  Bradley ,  "Stability  of  Subhawionic 
Limit  Motions  of  a  Slightly  Asymmetric  Missile , "  Ballistic  Research 
Laboratories  Memorandum  Report  No,  2494,  June  1975,  AD  B005079 . 

34.  C.H.  Murphy  and  J.W.  Bradley,  "Nonlinear  Limit  Motions  of  a 
Slightly  Asymmetric  Re-Entry  Vehicle,"  '  TAA  Journal  13,  July  1975, 
pp.  851-857.  (See  also  Ballistic  Rest  eh  Laboratories  Report  No. 
1755,  January  1975,  AD  A005165.) 

35.  B.G .  Karpov  and  J.W.  Bradley,  "A  Study  of  Causes  of  Short  Ranges 
of  the  8"  T31?  Shell,"  Ballistic  Research  Laboratories  Report  No. 
1049,  May  1958,  AD  377548. 

38,  W.  Hathaway,  R.  Whyte,  K.K.  Cobb  and  R.A.  Fierro,  "Analysis  of 

Aerohallistia  Range  Data  of  20mn  Projectiles  with  a  Moving  Internal 
Component ,  "  Proceedings  of  the  Fourth  International  Symposium  on 
Ballistics ,  17-19  October  1978 ,  American  Defense  Preparedness 
Association. 
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These  phenomena  have  been  explained37-38  by  assuming  either  (1)  a 
forced  circular  motion  of  an  internal  part  about  the  axis  of  the  pro¬ 
jectile  (hula-hoop  motion)  or  (2)  a  forced  precession  of  the  spin  axis 
of  the  internal  part  about  the  spin  axis  of  the  projectile.  In  both 
cases,  only  the  Fourier  component  of  the  motion  at  the  higher  coning 
frequency  of  the  projectile  was  considered.  Thus,  a  resonance  was  con¬ 
sidered  for  which  the  amplitude  of  the  internal  motion  was  constant. 

The  theory  predicted  a  relationship  between  the  growth  of  the 
higher  frequency  mode  and  the  spin  decay  and  in  both  cases  excellent 
quantitative  agreement  was  obtained.  Therefore,  the  derived  expres¬ 
sions  can  be  used  to  set  tolerances  for  manufacture  of  these  projec¬ 
tiles  . 


X.  ALMOST  SYMMETRIC  MISSILES 

Throughout  this  report  we  have  assumed  that  C  in  F:q.  (3)  was 

a 

zero  and  neglected  the  essential  asymmetry  present  when  the  pitch  and 
yaw  frequencies  are  not  equal  for  a  nonspinning  missile.  The  influ¬ 
ence  of  this  airplane-like  asymmetry  on  stability  of  a  spinning 
vehicle  was  studied  by  Phillips39  in  1948  and  developed  further  in  two 
recent  very  elegant  papers  by  Hodapp40-41.  The  concept  of  an  almost 
symmetric  missile42  was  introduced  in  1978  to  show  how  a  "little  bit" 


37.  W.G.  Soper,  "Projectile  Instability  Produced  by  Internal  Friction:," 
AIAA  Journal  16,  January  1978 ,  pp.  8-11. 

38.  C,H.  Murphy,  "Influence  of  Moving  Internal  Parts  on  Angular  Motion 
of  Spinning  Projectiles , "  Journal  of  Guidance  and  Con tro l  P, 

March- April  1978 ,  pp.  117-122.  ( See  also  Ballistic  Research  Lab¬ 

oratory  Memorandum  Report  No.  2731,  February  1977,  AD  A037338.) 

39.  W.H.  Phillips,  " Effect  of  Steady  Bolling  on  Longitudinal  and 
Directional  Stability,"  NACA  TN  1627 ,  June  1948. 

40.  A.E.  Hodapp,  Jr.,  "Effects  of  Unsymmetrical  Stability  Derivative 
Characteristics  on  Re-Entry  Vehicle  Trim  Angle  Behavior,  " 

Journal  of  Spacecraft  and  Rockets  11,  May  1974,  pp.  300-307 . 

41.  A.E.  Hodapp,  Jr.,  "Effects  of  Unsymmetrical  Stability  Derivative 
Characteristics  on  Re-Entry  Vehicle  Transient  Angular  Motion,  " 
Journal  of  Spacecraft  and  Rockets  13,  February  1976,  pp.  82-90 , 

42.  C.H.  Murphy,  "Angular  Motion  cf  Spinning  Almo st-Symmetric  Missiles,  " 
Journal  of  Guidance  mid  Control  2,  November-December  l'd?'J ,  pp. 
504-510.  (See  also  USA  ARRADCOM  Ballistic  Research  Laboratory 
Technical  Report  No.  02121,  November  1978 ,  AD  A063538.) 
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of  this  essential  asymmetry  would  affect  the  classical  tricyclic 
theory  of  a  symmetric  missile  with  aerodynamic  trim. 

If  is  not  zero,  the  pitch  frequency  of  a  nonspinning  missile, 
a 

is  not  equal  to  its  yaw  frequency,  •  If  the  spin  rate  is  out¬ 
side  the  region  between  these  two  frequencies  (the  resonance  region) , 
the  motion  can  be  described  by  a  pentacycle: 


i4>-  i(f4 

i  =  k  e  1+ke  ^  £  k  e  +ke  +ke 

1  2  3  4  5 


(29) 


where  <j/  =  2<f/  -  <f/  ;  ^  =  2<j/  - 


Mien  the  spin  is  in  the  resonance  region. 


s  (klR 


As  ,  ,  -As  _  . 

'  *  “<,R  "  '  h  \ 


u  .  t  „i*2  .  n 


w* . 


k  e  r2  ♦  v  e  5  (30) 

2  E 


According  to  Eq.  (30),  the  motion  grows  exponentially  when  the  spin  is 
in  the  resonance  region.  Another  unpleasant  characteristic  is  the 
existence  of  peak  values  of  |k  |  at  the  endpoints  of  the  region 


( ♦'  =  C  *'  a 


Four  specific  characteristics  of  the  motion  of  almost  symmetric 
missiles  [  | |<<|CM  I  ]  are : 
a  a 

(1)  The  genera]  motion  is  well  approximated  by  a  symmetric  missile 
with  average  coefficients. 
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(2)  Far  from  zero  spin  or  resonance  spin  rates,  the  first 
observable  modification  of  the  usual  tricyclic  motion  for  an  almost- 
synunetric  missile  is  the  appearance  of  a  2<f>/  -  <j/  frequency,  followed 

by  the  appearance  of  a  2 <p:  -  tp'0  frequency  as  the  asymmetry  becomes 
greater. 

(3)  Near  zero  spin,  both  cf  these  additional  frequencies  have 
substantial  amplitudes,  and  near  resonance,  the  2<j/  -  frequency 

has  a  substantial  amplitude. 

(4)  For  spin  in  the  resonance  region,  large  trims  and  exponential 
undamping  are  possible. 
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Figure  2.  Amplitude  plane  for  an  in-plane  cubic  damping  moment 

coefficient  three  times  larger  than  the  out-of-plane  cubic 
damping  moment  coefficient  (a  =  ?,  HQ  <  0  in  Eqs. 

(10-12)).  The  circular  limit  motions  are  stable,  but  the 
planar  limit  motion  is  unstable. 


the  gyroscopic 


Figure  5.  Amplitude  plane  for  a  cubic  Magnus  moment  (xj0  a2Q  >  0  in 

Eqs .  (20-21)).  The  separatrix  (dashed  curve)  divides  all 
initial  conditions  into  two  families:  those  that  yield  a 
motion  tending  to  zero  amplitude  and  those  that  yield  a 
motion  that  grows  without  bound. 
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savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


